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1. Introduction 

The two Higgs doublets model (2HDM) is the minimal extension of Standard Model (SM). 
The introduction of extra Higgs doublet can lead to the tree level flavour-changing-neutral-currents 
(FCNC), which should be strongly suppressed. Therefor it's necessary to restrict the general two 
Higgs doublets model in the Yukawa couplings sector. One of common ways is to impose Z2 
symmetry for the two Higgs doublets on the Lagrangian [jl]]. Under this symmetry there is only one 
Higgs doublet is allowed to interact with each type of fermions so that Yukawa coupling matrix is 
diagonal at any energy scale. 

Recently a more flexible way of avoiding tree level FCNC has been proposed [^. This idea 
simply requires that the two Yukawa couplings to the Higgs doublets should proportional to each 
other, so they can be diagonalized simultaneously. This restriction is fine in giving energy scale but 
the tree level FCNC will be reintroduced in higher energy scale ||3|]. There are also more general 
proposals which the tree level FCNC Yukawa couplings are suppressed enough, e.g. the Cheng- 
Sher ansatz [g]. 

We will study the properties of all these models by taking into account the theoretical and 
experimental constraints on FCNC. And then we will evolving all the Yukawa couplings to higher 
energy according to the Renormalization Group Equations (RGE). In some cases the off-diagonal 
matrix elements which relate to the FCNC grow quickly that indicates the certain assumptions 
behind the theory are not stale. Those theories are either fine-tuned or incomplete in certain way, 
e.g. there may be additional particles appearing when going to high energy scale. 

2. The general 2HDM and RGE equations 

One of standard convention to write the two Higgs doublets with the Goldstone bosons is 



Where and G° are the Goldstone bosons to be eaten by the EW gauge bosons during EW 
symmetry breaking, and is the chai^ged Higgs boson. The neutral Higgs scalar- can be divided 
into CP even scalars {h,H) and CP odd pseudo-scalar A. a and j8 is the mixing angle between 
Qi, H) and the two vacuum expectation values (VEV). 

For Yukawa coupling analysis, it's convenient to use Higgs basis 



Where the nonvanishing VEV is only assigned to Hi which plays the role of Standard Model Higgs 
doublet while H2 contains the new partciels and A. The general Yukawa coupling is 




(2.1) 



(2.2) 



Hi = cosj8<I>i+sinj3e"'*^a>2, 
H2 = -sinjSOi +cosj3e"'^02. 



(2.3) 



= QlHikI^Ur + Qj^HiK^Dr + LlHi^Er 

+QlH2PoUr + QlHipEDr + LlH2P^Er + h.c. . 



(2.4) 
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Where the subscribe "0" stands for flavor basis. In mass basis, the 3x3 matrix is related to 
diagonal fermion mass matrix by using bi-diagonalizing with the unitary matrices V[ , 

= V14VP = ^.£! (2.5) 

The coupUng matrix is still general and complex if there are no further restrictions. The Cheng- 
Sher ansatz suggests 

Pij - h • <-2.6) 

Where the m,- are the different fermion masses, the are expected to be of ^(1) and should be 
small enough to suppress FCNC to the observed level. In EW scale, the Z2 symmetric and aligned 
models can be teated as special case of Cheng-Sher ansatz. 



3. Numerical analysis 

3.1 SM input and Constraints 

The most stringent constrains for tree level FCNC is from neutral meson mixing. The master 
formula for tree level F° — F° mixing can be found in [|5|]. 

Using current experimental and theoretical data, we estimated the bounds for nondiagonal Xfj 
with the following assumption of Higgs scalar- mass: 

• mh=mH = mf, = 120 GeV: < 0.13 < 0.08 ,1,,^ < 0.03 , A,^ < 0.05}; 

• mh = mH = m GeV, ma = 400 GeV: < 0.30 , A^, < 0.20 , X^b < 0.08 , A,* < 0. 12 }. 

According to these results, we choose Xfj ^0.1 as a representive value which will be used as 
generic value later in RGB analysis. 

3.2 Z2 Symmetric Models 

The first example we study is the models with Z2 symmetry, in which the tan j8 is a physical 
parameter. The tree level FCNC is protected by the Z2 symmetry, so the nondiagonal Yukawa 
couplings stay as zero in any energy scale. The only thing we can study is to detect the place of 
Landau pole, where at least one of the Yukawas blow up. Beyond the Landau pole the perturbation 
theory is not valid anymore. 

In talbe |l| we show the diagonal X-j in terms of tan /3 for the four different 2HDM types in Z2 
symmetric models. The position of the Landau pole depends on the initial value of tan /3 in EW 
scale. In Fig. |l] we plot the upper and lower limit of tan /3 . The plots can be understood by whether 
the evolution is driven by Xn, Xhb -Xn or combination of them. 
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Type 




A,? 




I 


I/tanj8 


l/tanj3 


l/tanjS 


II 


l/tanjS 


-taiijS 


-tanjS 


miY 


l/tan/3 


-tan/3 


l/tanjS 


IV/X 


l/tan/3 


l/tanj3 


-tanjS 



Table 1: The diagonal X[- in 2HDM models with Z2 symmetry. 

Lower limits of tanp Upper limits of tanp 

Figure 1: The starting 
value of tan /3 as a function 
of the position of the cor- g.s 
responding Landau pole 
(ALandau-poie) in the differ- ,^^(3 
ent 2HDM types. There 

- 0.4 

are lower limits on tanp 
for all four types (left), but 0.2 
only type II, type III/X and 
type IV/Y have an upper ° 
Umit of tan )3 (right). Vandau-poie [GeV] ALan^jg^i-poie [GeV] 




3.3 Z2-breaking Models 

3.3.1 Aligned models 

In Yukawa alignment model, the two Yukawa couplings and are proportional to each 
other so they can be diagonalized simultaneously. In this model the is also diagonal in EW 
scale. However the nondiagonl element will start to grow when the Yukawas evolve to higher 
energy via RGE. Similarly to the Z2-breaking models, we limit the alignment model with three 
different cases: 

• Aligned I/II:A;f, A,f = 

• Aligned III: , Xfj = 

• Aligned IV: , = 

In Fig. ^ we plot the energy scale at which the Landau pole or large nondiagonal X^^j = 0. 1 is 
encountered as a function of pairwise combinations of the starting values for XV and A,f . 

3.3.2 Diagonal models 

Next we consider the models with Z2 symmetry breaking in either the up or the down sector. 
First we break the Z2 symmetry in the up-sector with X^ = X^ (type I) or X^ = — l/Xtt (type II), 
and set Xuu = Xcc- Then we break the Z2 symmetry in the down-sector with Xbb = X-/ (type I) or 
Xbb = -^/Xjf (type II), and set Xdd = A,*- 

3.3.3 Non-diagonal models 

In the end we consider the models of Z2 symmetry breaking from having non-zero non- 
diagonal elements in the up or down sectors. We set X^i- = 0.1 or XR: = 0.1 at the EW scale 
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A,°= A,'- or A,^= A,'- 



Figure 2: The constraints from the Landau pole 
and non-diagonal ^[^j- The plot shows the same 2 
results for Aligned I/II or Aligned III. The Xf( 
only play a very minor role so we don't show 
the plots for Aligned IV case. The areas inside aA. 
given contour are allowed by the requirement of 
the two condition above. The different contours ^ 
are as follows starting from the center; 10 
I0'°, 10^ 10^ and 300 GeV. The Z2 symmetric 
case is also plotted as a reference. -4 
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Figure 3: The energy scale where the Landau pole or non-diagonal Xf^j = 0.1 is reached as a function of 
Xss and Xht,. In the left (right) panels X^ = Xi,b{—1 /Xht,) and in all cases X^ = Xhb- 
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Figure 4: The constraints from the Landau pole and the off-diagonal elements as a function of Xjf and 
the off-diagonal elements Xj^j (up) or Xj^^j (down) at the input scale for the type I (left) and type II (right) 
relations for the diagonal elements. 

and then evolve to higher energy. The analysis shows the constraints from off-diagonal elements 
are small compare to the previous cases. 

4. Conclusions 

RGE evolution is a useful tool to analysis different 2HDM models on stability of underlying 
assumptions. A quickly appearance of Landau pole or large off-diagonal Yukawa coupling un- 
der RGE evolution may indicate the model is either fine tuned or incomplete, e.g., new particles 
appearing in high energy. 
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